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Abstract. 

The orbital dynamics of Cyg X-3 are a key to understanding this enigmatic X-ray 
binary. Recent observations by the RXTE ASM and the OSSE instrument on GRO 
enable us to extend the baseline of arrival time measurements and test earlier models 
of orbital period evolution. We derive new quadratic and cubic ephemerides from the 
soft X-ray data (including ASM) . We find a significant shift between the predicted soft 
X-ray phase and the light curve phase measured by OSSE from ~ 44 to 130 keV. Some 
of the apparent phase shift may be caused by a difference in light curve shape. 



INTRODUCTION 



Cygnus X-3 is a unique and poorly understood X-ray source which may represent 
a very short-lived, transitional state of binary evolution. A study of its dynamics 
(via the orbital ephemeris), and local environment (via the light curve and spec- 
trum) may have broad implications for X-ray binaries in general. The period (4.8 
hr, presumably orbital) is characteristic of low mass X-ray binaries. However, there 
is evidence that the companion is actually a high- mass Wolf-Rayet star [1]. The 
mass loss from the stellar wind of such a star could explain the large positive P 
measured in Cyg X-3, which is inconsistent with mass transfer via Roche lobe 
overflow [2] . The light curve, observed in IR and X-rays, is asymmetric with a non- 
zero eclipse, resulting, in most models, from interactions in a dense cocoon or wind 
around the system [3-6]. If the asymmetry is instead caused by an elliptical orbit, 
the light curve shape is expected to change significantly over time due to apsidal 
motion [7] which could also produce some or all of the apparent period increase [8] . 
Therefore, measurements of the light curve shape are directly related to the orbital 
dynamics. 



OBSERVATION SUMMARY 



OSSE [9] observed Cyg X-3 for ~ 65 days in 7 intervals between 1991 May 30 and 
1994 July 12. We use the 2-min background-subtracted ~ 50 keV to 10 MeV spectra 
produced by the standard OSSE spectral analysis routines for these 7 observations. 
There was no significant flux contribution from Cyg X-l in any of the source or 
background pointings. Some results of the earlier OSSE observations have been 
described elsewhere [10]. 

For the ASM we used the background subtracted 2-10 keV Cyg X-3 counting rates 
from the quick-look results provided by the ASM/RXTE team over the WWW. The 
data analyzed here cover 1996 Feb 22 - 1997 Apr 11. 

TABLE 1. ASM Arrival Times 

Interval Cycle Arr. Time ot 

(HJD - 2440000.) (HJD - 2440000.) (days) 

10135.8- 10235.4 46251.0000 10185.6237 1.0 x 10~ 3 
10236.0-10335.9 46753.0000 10285.8687 1.1 x 10~ 3 

10335.9- 10435.8 47254.0000 10385.9134 1.4 x 10~ 3 
10435.8-10535.7 47754.0000 10485.7588 1.2 x 10~ 3 



ASM ARRIVAL TIME DERIVATION 

To correct for large, long term changes in the source DC flux level the ASM 
2-10 keV Cyg X-3 rates were first "cleaned" by fitting a line to successive 10-day 
segments of the data and subtracting the fit. The linear fit removes trends on time 
scales much longer than an orbital period. This process did not significantly change 
the arrival times derived from the data. 

The cleaned data were divided into approximately 100-day intervals. Long inter- 
vals are needed to average over the cycle-to-cycle light curve (LC) variations and 
to completely sample the Cyg X-3 orbit. The data in each interval were then cor- 
rected to the SSB and epoch-folded (using a constant period ephemeris) to produce 
orbital light curves. We fit the resulting light curves with the EXOSAT template 
[11] to determine the relative phase of the LC minimum. This phase is then used 
to calculate an arrival time for a cycle near the middle of the observation interval 
[12]. The arrival times are shown in Table 1 in heliocentric JD (HJD) — 2440000. 

Figure 1 shows differences (residuals) between the arrival times predicted by a 
constant period ephemeris and those observed by the ASM and by other soft X-ray 
instruments going back to 1970 ( [13] and references therein, plus [14]). Overplotted 
are fits with quadratic and cubic ephemeris models. 



0.06 



0.04 



3 0.02 



CO 

>. 

CO 

;o 
_co 

CO 



■o 

CO 
CD 

DC 



0.00 



-0.02 



oASM 
xOSSE 

^Ginga, ASCA [13] 



Const ephem 
Quad, excl. OSSE 
Cubic, excl. OSSE 
Cubic (ref. 11) 




2.0*1 0° 4.0*1 0° 6.0*10° 8.0*10" 
Time (JD - 2440000) 



1.0*10 



1.2*10 



FIGURE 1. Arrival time residuals for OSSE, ASM, and historical data. The errors on the ASM 
points are approximately equal to the symbol size. Also shown are the best-fit quadratic and 
cubic ephemerides. An extrapolation of an older cubic fit [11] with a larger negative P is shown 
for comparison. 



ORBITAL EPHEMERIS 



The ASM arrival times were combined with those from earlier soft X-ray obser- 
vations and fit with both quadratic and cubic ephemeris models (Figure 1). The 
resulting best-fit parameters are shown in Table 2. Neither quadratic nor cubic 
models provide a statistically acceptable fit to all the data. Possibly the errors 
in earlier data have been underestimated due to cycle-to-cycle variations or other 
systematic effects [13]. No single data set drives the bad x 2 . 

The cubic term is still statistically significant, but the fitted magnitude has de- 
creased with time: from P = 1.2 x 1CT 10 [11] to 4.1 x 1CT 11 yr 1 in this work. 
The wind model predicts essentially zero P. Future ASM data should provide a 
definitive test of the reality of the non-zero P. 



OSSE ARRIVAL TIME DERIVATION 

Following the same procedure used for the ASM data, arrival times were de- 
termined for the each of the 7 OSSE observations of Cyg X-3 from 1991 to 1994 
using light curves of 2-min background-subtracted counting rates ~ 44-130 keV. 



TABLE 2. PRELIMINARY Best-Fit Quadratic and Cubic Ephemerides 



T n = T Q + P n + c n 2 



T n = T + P Q n + c n 2 + dn 3 



To 
Po 



2440949.89185 ± 0.00088 HJD 
0.199684393 ± (8.9 x 10~ 8 ) d 
(6.06 ±0.17) x 10- 11 d 



2440949.8967 ± 0.0014 HJD 
0.199683305 ± (2.6 x 10~ 7 ) d 
(1.17 ±0.13) x 10~ 10 d 
(-7.54 ±1.71) x 10~ 16 d 
136.37 (83) 
2.0 x 10~ 4 



Co 
d 



X 2 (dof) 



168.42 (84) 
1.3 x 10~ 7 



X 2 (dof) 



Conf. level 



Conf. level 



The arrival times for the 7 observations were compared with the values predicted 
using the soft X-ray ephemerides (Table 2). The OSSE points fall systematically 
below the predicted curves on the residual plot. 

The OSSE phase minimum (~ 44-130 keV) is significantly (6- 7a) earlier than the 
predicted soft X-ray (1-10 keV) minimum for both quadratic and cubic ephemerides 
(At ~ — 20±3 min, A0 ~ — 0.065±0.010). The data are consistent with a constant 
offset over more than 3 years of observations. 



Light curves (Figure 2) were produced from the ASM and OSSE data using the 
quadratic ephemeris in Table 2. The plotted errors in the ASM data reflect only 
the reported statistical errors. The light curves were fit with the EXOSAT soft 
X-ray template [11], varying phase, amplitude, and DC intensity to determine the 
phase of minimum and the consistency of the overall shape with the template. 

Both the OSSE and ASM light curves are statistically inconsistent with the 
canonical X-ray template [11]. Qualitatively, however, the template reasonably 
describes the ASM data. The large y 2 is due to the fact that the actual observed 
variations in each bin are larger that the statistical errors due to underlying source 
fluctuations. 

Evolution of the light curve shape with time is expected if the apparent orbital 
period change is due to apsidal motion [7,8]. There is no evidence for this in the 
ASM observations. 

The OSSE LC appears somewhat more symmetric than the template, with a 
faster rise, but this is difficult to constrain with current data. The differences in 
LC shape may contribute to the apparent phase shift with energy. 



ORBITAL LIGHT CURVE ANALYSIS 
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FIGURE 2. ASM (top) and OSSE (bottom) light curves with best-fit X-ray templates over- 
plotted (dashed line). The dashcd-dotted line in the OSSE plot shows the best fit to the OSSE 
data with the ASM phase. The shift between the two energies is apparent. 
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